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Magnetars!, a population of isolated neutron stars with ultra- 
strong magnetic fields of ~ 10'* — 10'° G, have been increasingly 
accepted to explain a variety of astrophysical transients. A nascent 
millisecond-period magnetar can release its spin-down energy and 
power bright sources such as Gamma-ray Bursts? (GRBs) and 
their subsequent X-ray plateaus*, Super Luminous Supernovae 
(SLSNe)*>, and the fast X-ray transients such as CDF-S XT-2°. 
Magnetars with ages of 10° — 10* years have been observed within 
the Milky Way Galaxy, which are found to power diverse transients 
with the expense of their magnetic energy, in the form of giant 
flares and repeated soft-y-ray or hard X-ray bursts”* and occa- 
sionally fast radio bursts (FRBs)*°. Magnetar giant flares were 
also detected as disguised short GRBs from nearby galaxies 1128, 
Here we report the identification of a GRB as a hyper flare of mag- 
netar in a nearby galaxy. The magnitude of the hyper flare is about 
one thousand times brighter than that of a typical magnetar giant 
flare. A significant ~ 80 millisecond period is detected in the de- 
caying light curve. Interpreting this period as the rotation period 
and given a magnetic field strength typical for a young magnetar, 
the age of the magnetar is constrained to be only a few weeks. The 
non-detection of a (superluminous) supernova nor a GRB weeks 
before the event further constrains that the magnetar is likely born 
from an off-axis merger event of two neutron stars. Our finding 
bridges the gap between the hypothetical millisecond magnetars 
and the observed Galactic magnetars, and points toward a broader 
channel of magnetar-powered gamma-ray transients. 

Recent observations of the extragalactic magnetar giant flare 
(MGF)!!!2!3, GRB 200415A, suggested that MGF GRBs follow the 
third track (other than long and short GRBs) in the rest-frame peak 
energy vs. isotropic energy (Ep — Eiso) plane!?. Motivated by this 
new track, we performed a systematic search in the Fermi/GBM GRB 
archives '*, aiming at finding additional MGF GRB cases. Our search 
started from looking for those GRBs with relatively high E, and low 
Eiso (or low fluence if the redshift is not available), so they can be out- 
liers of the short and long GRB populations in the E, — Eiso diagram. 
Our search quickly returned a strong candidate, GRB 130310A (Figure 
1). The burst consists of two main emission episodes: a precursor and a 
main burst, lasting a total of ~ 4.2 seconds (see Methods and Extended 
Figure 1). The precursor triggered Fermi/GBM © at 20:09:41.503 on 
March 10” 2013 UTC (hereafter To) and lasted for approximately 0.8 
s. It is characterized by a thermal spectrum with kT = 45.06*337 
keV (Methods). The main burst occurred at about To + 3.8 s and pre- 
sented a sharp peak followed by a series of erratic overlapping pulses 
and lasted for approximately 1.3 s (Methods). Such a profile broadly 


resembles the previous observations of magnetar giant flares !6, which 
are characterized by a hard, sharp spike followed by a soft, long-lasting 
tail. The time-integrated spectrum of the main burst is characterized by 
a Cutoff Power Law (CPL) model with E, = 273260 "er a keV and 
a= —1.10400. Such an E, value is already significantly higher than 
most GRBs, making GRB 130310A, even with any assumed redshift 
between 0.0001 and 10, a distinct outlier from the long collapsar-type 
GRB track in the Ep — Eiso diagram (Figure 2). In consideration of 
its strong spectral evolution (Methods) as well as its spike+tail feature, 
we treat the 27-ms sharp peak (SP) (ranging from 4.108 to 4.135 s; as 
indicated in Extended Data Figure 1) as the characteristic emission of 
the event, and the following emission as its extended radiation. The 
time-integrated spectrum of the SP is best fitted by a CPL model pa- 
rameterized by œ = —1.19 + 0.02 and E, = 9.8122 MeV, with a fluence 
of 9.88+020 x 10-76 erg cm™ (Methods). This suggests that the bulk 
Lorentz factor should be at least ~ 430 (Methods), being the highest 
among all observed MGFs observed so far!!. 

Since there is no reportedly redshift measurement of GRB 
130310A, we assigned z as a free parameter ranging from 10~* to 10, 
and overplotted the corresponding values of the rest-frame E, and Eiso 
of the SP onto the Ep-Eiso diagram, as shown in Figure 2. Interestingly, 
due to its high E,, the event lies far away from typical long and short 
GRBs tracks. Instead, it is consistent with the third track for the MGF- 
GRB population. The redshift of the GRB is constrained in the range 
of [0.0079, 0.0595], assuming that the burst follows the MGF track. 
We further searched for the host galaxy (Methods) of the burst in the 
6df Galaxy Survey (6dfgs) database !™!8 within such a redshift range 
inside the overlapped region between the LAT error circle (Methods) 
and the IPN error box!®. Our search yields only one galaxy within 
that redshift range and error box, which is g0927191-170053 located 
at RA = +09'27"19.70° and DEC = —17°00'53.0” with a redshift of 
0.0155. At such a redshift, the isotropic energy (peak luminosity) of SP 
is 5.61*95' x 10°8 erg (2.307952 x 10°° erg s™!). With such a low energy 
and a high spectral peak, GRB 130310A consistently lies at the high- 
end of the MGF GRB population track (Figure 2). The total isotropic 
energy of the main burst is 2.377918 x 10 erg, which is around a thou- 
sand times higher than that of GRB 200415A. We hence consider it as 
a “hyper” flare of a magnetar in this study. 

The magnetar nature of GRB 130310A is further manifested by the 
period detection in the light curve. We employed the Lomb-Scargle 
(LS) method?°?! to search for periodic signals in different time win- 
dows with different time scales in different energy ranges (Methods). 
As shown in Figure 1, we performed the LS calculation on the original 
light curves as well as the detrended light curves resulted from three 
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Figure 1 | Light curve and periodic signal detection of the hyper flare event, GRB 130310A. 
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middle blue and lower right black solid lines show the best-fit correlations for the 
MGF, short and long GRB populations, respectively. Dashed borderlines show the 
30 regions for each correlation. The red dashed line represents the locations of 
the sharp peak of GRB 130310A at different redshift values. The red square marks 
its location at a redshift of z = 0.0155, which is determined by its candidate host 
galaxy, g0927191-170053. The two red circles mark the allowed edge if the MGF 
folows the MGF-GRB track. All error bars represent 10 uncertainties. 


different detrending methods. Our analysis yielded a significant peri- 
odic signal of 12.4 Hz with a confidence > 5o found in the 8-1000 keV 
light curve in the time region between 4.4 s and 4.9 s (Methods). Such 
a period is unprecedented in the magnetar observations and has never 
been observed in GRB events. 

The period of ~ 80 ms can be directly connected to the neutron 
star rotation. Assuming magnetic dipole radiation, the spin-down rate, 
P, can be estimated from the spin period, P, and its estimated surface 
magnetic field strength, B,, as”: 
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where c is the speed of light, Mo is the solar mass, 7 ~ 2MR2/5 is 
the moment of inertia, M and R, are the mass and surface radius of 
the magnetar, respectively. We normalize M to be 2Mo considering 
that the magnetar is likely born from a binary neutron star merger, as 
explained below. Here we have adopted Q = 10”Q, in cgs units. The 
magnetic field of a young megnetar can be as high as~ 10!° G? and 
is normalized to 10!55 (or 3.2x10!5) G in this study. The characteristic 
age of the magnetrar can be estimated as”: 
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which is about three weeks for nominal parameters. For simplicity, 
we further assume that there is about 50% uncertainty in assigning the 
B, value, so the minimal and maximum age of the magnetar can be 
calculated as Temin = 11 days and Te max = 56 days respectively. 


(2) 


At age T,, the spin-down luminosity is 
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We note that such an Lsa is much smaller than the observed peak lumi- 
nosity of GRB 130310A, suggesting that the hyper flare is not powered 
by the spindown of the magnetar, but is rather powered by significant 
magnetic energy dissipation in this early phase of a magnetar’s life. 
One may compare the GRB energy with the total energy available 
from the source. Considering a beaming factor of f ~ 0.01, the total 
energy of GRB 130310A may be estimated as Ey = f Eiso = 2.377018 x 


-0.18 
10” f> erg. The total spin energy of the magnetar at the epoch is 
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and the total magnetic energy of the magnetar can be estimated as 
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One can see that the GRB energy is within the energy budget of the 
magnetar if f < 1. The hyper flare likely originates from the instanta- 
neous release of the magnetic energy due to the dissipation of magnetic 
energy. The emission is likely highly beamed so that magnetar rotation 
can leave a periodic signal in the light curve. 

The estimated position of the magnetar powering GRB 130310A 
in the P — P diagram is presented in Figure 3. One can see that it 
bridges the hypothetical millisecond magnetars born in GRBs/SLSNe 
and Galactic magnetars. 

We searched for potential GRBs from the archival data as the pos- 
sible progenitor of the magnetar. Our search window covers the time 
range between To — Tenax and To — Temin: We followed the burst search 
method described in Ref.”>. The search yielded no significant signal 
in either triggered or untriggered GRBs samples around the location of 
GRB 130310A. The duty cycle, calculated by counting the time-span 
during which the GRB 130310A was in a good field of view (FOV) of 
Fermi/GBM,, is 64.6%, suggesting that we had about 35.4% probability 
to miss the progenitor GRB event if there was any. Assuming a typical 
GRB spectrum parameterized by a CPL model with a = —1 and E, 
= 300 keV, the 10-s flux upper limit in 10-1000 keV at To-23.9 day is 
7.01 x 1078 erg cm~? s~!, corresponding to a luminosity upper limit of 
3.98 x 10% erg s~!, which is ~ 6000 times fainter than the peak lumi- 
nosity of GRB 130310A. Our search suggests that a progenitor GRB 
must be significantly off-axis if it was not missed during the 35.4% off 
time. 

We further utilize the archival data in the optical band to constrain 
the progenitor type of the magnetar. No supernova discovery was re- 
ported from the host galaxy within To+ 6 months. The Panoramic Sur- 
vey Telescope and Rapid Response System (Pan-STARRS)”® largely 
covered the region of the magnetar from 2012-06 to 2014-07. We 
searched the PAN-STARRS data archive’ and found that there were 
15 observations, from November 04, 2012 to November 18, 2013, with 
a total exposure time of ~ 912 seconds, covering the host galaxy re- 
gion, which can place a series of 30 non-detection upper limit for a 
point source around To (Methods), as plotted in Figure 4. Such upper 
limits can rule out the existence of a Type-Ic supernova associated with 
a long GRB or an SLSN at any time later than To — Te,max days (Meth- 
ods). On the other hand, a kilonova or an off-axis afterglow of a GRB 
170817A-like event at ~ 70 Mpc, can survive these constraints. Even 
for an on-beam short GRB, the afterglow emission may also avoid de- 
tection if it occurred in the gaps between PAN-STARRS observations. 
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Figure 3| P-P diagram. All known neutron stars collected in different categories 
in the ATNF Pulsar Catalogue?” are plotted with different symbols as indicated in 
the top-right corner. The weeks-old magnetar in this study is overplotted with 
a red star. The hypothetical newborn millisecond magnetars with B~ 10!5-!6 G 
are located in the diamond region. The arrow indicates the evolution path from 
a newborn magentar to a weeks-old magnetar, assuming that B does not decay 
significantly. 
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Figure 4 | Optical upper limits from the archival data and their constraints 
on optical emissions from the magnetar progenitor. The arrows show the multi- 
band 30 upper limits for the non-detection of a point source transient in the skirt 
of the galaxy g0927191-170053 constrained from the PANN-STARR archival data 
(Methods). The solid green line shows a typical superluminous supernova light 
curve, assuming that it occurred at To - 23.9 day at ~ 70 Mpc. The brown solid 
line shows a typical type Ic supernova light curve assuming that it occurred at 
To - 23.9 day at ~ 70 Mpc. The orange, purple and navy curves stand for the 
kilonova, on-axis and off-axis optical afterglow light curves, respectively, for a 
GRB 170817A-like event located at 70 Mpc. 


All these suggest that the progenitor of the magnetar is likely a neutron 
star merger event, with the off-axis viewing angle preferred because 
of its much higher probability. We note that the merger origin is also 
consistent with the fact that its host, g0927191-170053, is an old-type 
elliptical galaxy 8. 

To date, the age of the magnetar is approximately 8.5 years old. 
Using the total energy of 2.37 x 10” erg derived above and the stan- 
dard afterglow parameters (Methods), we calculate the afterglow flux 
as ~ 1.02 x 107 / 3.04 x 10°° mJy in the X-ray/optical bands (Meth- 
ods), which are way below the detection thresholds of the current mis- 
sions/telescopes in operation. One possible way of detecting emission 
from the system may be fast radio bursts if the magnetar can survive 
until now. At ~ 70 Mpc, an FRB with a typical luminosity of 10” erg 
s7! should be readily detectable with most radio telescopes. 


METHODS 


Fermi data analyses 

The hyper flare event, GRB 130310A, was a bright burst detected 
by several high-energy missions including Fermi/GBM!>, Fermi/LAT”?, 
Konus-Wind*°, and Suzaku WAM?!. We performed our analysis based 
on the time-tagged events (TTE) data collected by the Gamma-ray Monitor 
(GBM) onboard the Fermi satellite. The data were obtained from the Fermi 
Public Data Archive (https: //heasarc.gsfc.nasa.gov/FTP/fermi/ 
data/). We use the data of sodium iodide (Nal) detectors n9 & na and 
germanium oxide (BGO) detector b1, which have the smallest angular sepa- 
rations with respect to the location of GRB 130310A and the highest bright- 
ness. We also utilized the Large Area Telescope (LAT) data covering 30 
MeV to 100 MeV (LLE) and 100 MeV—300 GeV energy ranges. The stan- 
dard data reduction procedure was described in3?3334, The key results of 
our analysis are outlined below. 

Light curves and duration: Light curves in different energy ranges are 
extracted from the TTE data following the procedures described in Ref**. 
Those light curves are searched for possible periodic signals. The multi- 


wavelength light curves with typical energy ranges are presented in Extended 
Data Figure 2. The burst duration (Too) is calculated in the standard range of 
8 keV- 1 MeV. As shown in Extended Data Figure 1 and listed in Extended 
Data Table 1, the To9, defined by the time interval during which 90% of total 
counts are detected, is ~ 1.3 s for the precursor and is ~ 2.93 s for the main 
burst. 

Time lag: Time lag measures the difference of photon arrival time in 
different energies *> and is often used as a probe to infer the emission size. 
We utilized the cross-correlation function (CCF)**”® to calculate the time 
lags of light curves among different energy bands between Ty + 4.1 s and 
To + 4.6s, following the method described in ref. 37. As shown in Extended 
Data Figure 2, no significant lag is detected for the light curves in different 
energy bands. The zero-lag result suggests that the emission region of GRB 
130310A is small, consistent with the speculation that it is an MGF GRB. 

Amplitude parameter: The amplitude parameter f = ae is defined as the 
ratio between the peak flux F, and background flux Fg (Ref.*8). A small f 
often suggests that the signal is more likely affected by the ’tip-of-iceberg” 
effect. The f parameter of the sharp peak of GRB 130310A is 11.22 + 1.32, 
which is the highest among all the GRBs (Extended Data Figure 4). 

Spectral analysis: For the main burst ranging from To + 3.85s to 
To + 5.25s, we performed both time-integrated and time-resolved spectral 
analyses using the data from the two Nal and one BGO detectors as men- 
tioned ahead. The time slices for time-resolved analyses are obtained ac- 
cording to the light curve brightness profile in such a way that the pho- 
ton count of each spectral channel is greater than twenty. All those inter- 
vals are listed in Extended Data Table 2. The spectral files, including the 
total observed count spectra, the background spectra, and the detector re- 
sponse matrices (DRMs), are extracted from the event files using the same 
method as described in Ref. °”. Four frequently used spectral models, namely 
Band function (Band), Black body (BB), simple power-law (PL), and cutoff 
power-law (CPL), are employed to fit the observed spectra. In particular, the 
CPL model is defined as*° 


N(E) = AE” exp[-(« + 2)E/Ep], (6) 


where a is the power-law photon index, E, is the peak energy of the fitted 
spectrum in units of keV, and A is the normalization factor. 

The time-integrated spectrum of the main burst is characterized by a 
CPL model with E, = 2454.34*12)78 keV and a = —1.05+%01, We noticed 
that such an E, value is already significantly higher than that of most typical 
GRBs. 

The time-integrated spectrum of the shart peak alone is characterized 
by a CPL model with E, = 9802.23+723488 keV and œ = —1.19 + 0.02. We 
noticed that such an E, value is already significantly higher than that of most 
typical GRBs. The total fluence of the first peak is 9.887000 x 107 erg cm™. 

By comparing the goodness of the fits in each time slice of the main 
burst, we found that the CPL model is the best model that adequately de- 
scribes the observed data with the lowest Bayesian Information Criteria 
(BIC)!, we thus only employ the CPL model to perform the time-resolved 
spectral fitting. 

Interestingly, the time-resolved spectral analysis suggests even more ex- 
treme values of Ep. With significant spectral evolution within the To inter- 
val, Ep reaches at ~ 10*}3 MeV around the peak region of the light curve. 
As shown in Extended Data Figure 3, we note that a and E, evolution shows 
an intensity tracking behavior as suggested by 434, 

High-energy photons: We extracted the LAT data in the energy range 
of 100 MeV—300 GeV from To to To + 800 s. The region of interest 
(ROI) is a square with the width of 10° and the center at the position 
RA, Dec = 142.34, -17.23 (J2000). Using the standard Fermitools*, and 
with a cut on zenith angle at 100° and the instrument response function 
set P8R3_SOURCE_V3, an unbinned likelihood point source analysis is 
performed. The GRB spectral model is assumed to be Powerlaw2 ? with 
Emin = 20 and Emax = 20000. Consequently, the spectral index is fitted to be 
—1.9+0.1. Finally, we obtained 20 photons with different probabilities orig- 


'BIC4! is defined as BIC = —2InL + klnn, where £ is the maximum likelihood, k is 
the number of parameters of the model, and N is the number of data points used in the fit. 

2https://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/source_ 
models. html 


inating from GRB 130310A. Those high-energy photon events are plotted in 
Extended Data Figure 2a & 2b. 


Amati relation 

The correlation between the GRB isotropic energy Ey;iso and the rest- 
frame peak energy Ep, = (1 + z)Ep is commonly refereed to as Amati rela- 
tion*ć, and can be written as log Epz = a + blog Eyjiso. Typically, long and 
short GRBs, as well as MGF GRBs follow different tracks in the diagram 39. 
We plotted the E,, — Ey is. diagram (Figure 2a) using GRB samples with 
known redshift from Ref.*° and Ref.*’. Utilizing the MCMC method, the 
optimal fitting parameters and 3c" uncertainties are constrained. We obtain 
a = —25.01*331 and b = 0.52*}0% for long GRBs, a = —19.92*7/* and 


b = 0.45+%15 for short GRBs, and a = -13.20*7 8 and b = 0.354-12 for 
MGF GRBs. The hyper flare event, GRB 130310A, is a significant outlier 
from both the short and long GRB tracks, but is consistent with being an 
MGF GRB. Assuming that it follows the MGF track, the redshift of GRB 


130310A can be constrained within the range from 0.0079 to 0.0595. 


Periodicity Measurement 

To maximize the signal detection, we combined all the photon events 
from the three detectors n9, na, and nb to search for periodic signals. The 
light curves in different energy bands are all binned with 0.0005s, which 
ensure sufficient resolution for the periodicity search. We searched possible 
periodic signals in the light curve of the main burst by using the Lomb- 
Scargle (LS) periodograms method. In addition to the original light curves, 
our search also considers the detrended light curves, with the latter aiming 
to remove the effect of global temporal evolution of the event. The following 
three detrending models are employed and fitted to the search phase of the 
light curve: 


1. Exponential model f(t) = ae” + c, indicated as yellow solid line in 
Figure 1. 


2. Whittaker Smooth*®“? (WS) model, indicated as blue solid line in Fig- 
ure 1. WS is a fast and wildly used smooth approach, which can well 
fit the trend of the light curve and retain the periodic signal by giving 
an appropriate smoothing parameter. 


3. A model using the multivariate adaptive regression splines algo- 
rithm*°, coded as py-earth model and indicated as the purple solid line 
in Figure 1. The multivariate adaptive regression splines algorithm is 
a non-parametric regression method that builds multiple linear regres- 
sion models on each different partition, which can efficiently catch the 
global temporal features of the light curve. 


The LS periodogram calculated based on the original light curve as well 
as the de-trended light curves by the above three methods are shown in the 
right column of Figure 1. A significant periodic signal of 12.4 Hz, corre- 
sponding to a period of 0.08s, is detected in all four situations. The strongest 
detection of the signal lies in a energy range of 50-300 KeV, and time range 
of 4.4-4.9 s. The confidence level of the detection in all cases exceeds 50, 
with a false alarm probability (FAP)*!>? of < 5.7 x 1077. 

The periodic signal as well as its detectable range are also verified by 
applying the weighted wavelet Z-transform) to the detrended data, as shown 
and Extended Data Figure 7. 

The lack of detection of the periodic signal after 4.9 s in the extended 
tail is not unexpected as the observed flux can become less modulated by 
the underlying periodic signal when it decreases exponentially and emerges 
to the background level. This can be illustrated in Extended Data Figure 8, 
where we generated a series of Fast-Rising-Exponential-Decaying (FRED) 
shape light curves with injected noise similar to that of GRB 130310A, as 
well as a periodic f = 12 Hz signal, S ,, characterized by different amplitude 
A, time scale r, which is formulated by 


Sp) = sin(27 ft) (7) 


exp[@ - 44/7] 


As shown in Extended Data Figure 8, our simulations indicate that the 
periodic signal can become undetectable after 4.9 s in the extended tail under 
certain configurations of the S ,(¢). 


Implication on the bulk Lorentz factor 

Following the derivation from Ref.>*, for a photon spectrum modelled 
as fe~**, in which f is the normalization factor in unit of ph cm~? s~! keV~! 
and a, is the inverse number of the photon index, œ, in correspondence with 
the definition in Eq. 6, the lower limit of the bulk Lorentz factor can be 
constrained by: 
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where Emax is the photon with the highest energy, ôT is the minimal vari- 
ability time scale constrained through the Bayesian block method, d is the 
distance of the burst, or is the Thomson cross-section and 7(a@,) is the 
correction factor which reflects the averaging effect of the pair production 
cross section. The form of 7(a,) varies in the literature*4>° and we adopt 
n(a,) = (3/8) + a)! as discussed in Ref.>>. 

Observationally, the minimal variability time scale of GRB 130310A is 
7.5 ms. The sharp peak whose spectrum is fitted as a CPL model parame- 
terized by f=98.57, a, = 1.19, and Emax ~ Ep = 9.8 MeV. By putting those 
numbers in Eq. 7, we can calculate the lower limit of the bulk Lorentz factor 
as Ymin = 430. 


Localization with Fermi/LAT and IPN 

By removing the source model of GRB 130310A from the fitted model, 
we construct the Test Statistics (TS) map by gitsmap, and we localize GRB 
130310A at RA = 142.52 and DEC = -17.16 with the maximum TS value. 
To obtain the confidence interval, we subtract the TSmax from the TS map 
to build a Localization Test Statistics (LTS) map in Extended Data Figure 
5. The white, orange, green, and black ellipses denote the 68.3%, 95.5%, 
99.7%, and 99.9% confidence areas, respectively. The LAT counterpart of 
the peak of the GBM light curve is located with RA = 144.50 and DEC = 
-19.02. 

The InterPlanetary Network (IPN) reports a constraint of the burst lo- 
cation from a group of spacecraft equipped with gamma-ray burst (GRB) 
detectors !°. The overlapping region between the IPN error box and the 99% 
confidence of the LAT box is regarded as the most probable location region 
of the hyper flare event and is used for the host galaxy search. 


Host galaxy search 

Within the redshift range constrained by the Amati relation and spa- 
tial range of the overlapping IPN & LAT error boxes, we performed a host 
galaxy search in the 6df Galaxy Survey (6dfgs) database !™!8, which records 
more than 100,000 galaxies in the nearby universe. Only one host galaxy, 
20927191-170053 with RA = +09"27"19.07° and DEC = —17°00’53.0” 
and redshift = 0.015492, was found. The galaxy is also cross-matched in 
the Two Micron All Sky Survey Extended Source (2MASS), as 09271905- 
170052856. The location of the galaxy is marked in Extended Data Figure 
5. 


Optical upper limits of the magnetar emission around To 

We search for all available archival optical sky-survey data in the hope 
of finding observations that can cover several months around To. Those 
observations can provide upper limits, if not detection, of a point source 
associated with the magnetar’s optical counterpart. 

Our investigation yields one set of archival data obtained from the 
Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). 
Around To, Pan-STARRS covered the host galaxy of the magnetar for a total 
period of about one year, from November 04, 2012 to November 18, 2013, 
with a total observational time of 912 seconds”®. Those observations are 
listed in Extended Data Table 3, and are available as single-epoch “warp” 
images in Pan-STARRS|1 data archive, which are astrometrically and photo- 
metrically calibrated*. No variable point source was found within a circular 
region of 3’ radius from the galaxy center in those observations. For each 


3https://outerspace.stsci.edu/display/PANSTARRS/PS1+Warp+images 


warp image, we calculate the 30 upper limit of a merger-type point source 
detection using the following procedure: 


1. For each pixel in the warp image, we calculate its A-B magnitude fol- 
lowing the guide as listed in Pan-STARRS document* and Ref>”>° as 
(10) 


where A is the magnitude zero point, C is the photon counts in each 
pixel, texp is the exposure time in units of second. 


mag = A—2.5 X logC + 2.5 X logtexp, 


2. The K half-light radius is Rso = 4.94” °°. Since the magnetar is likely 
a merger event, we assume it occurred at the outskirt of the galaxy 
with the normalized offset rog with a similar value as in the case of 
GRB 050724, namely ror = Ror = 0.68. This allows us to place the 


magnetar at a distance of Rof = 3 36” from the galaxy center. 


3. With an same inclination angle of the host galaxy, we selecte a region 
of an ellipse annulus with the outer semi-major axis of 4.30” and the 
inner semi-major axis of 3.19”, as shown in Extended Data Figure 6. 
We then calculate the magnitude of each pixel within the annulus and 
plot their probability density distribution, as shown in Extended Data 
Figure 6. The 30 value of such distribution is used to determine the 30 
upper limit of a point source at radius Rog. 


Through the above steps, we obtained the 30 upper limit of a point 
source at r for a total of 15 warp images, as listed in Extended Data Table 3 
and plotted in Figure 4. We note that 8 warp images in Extended Data Table 
3 are considered “bad” (e.g., the magnitude value of each pixel is marked as 
“nan” due to gaps between detectors, bad pixel regions, etc*.) so their upper 
limits are unavailable. 


Possible optical light curves of the progenitor 

The progenitor of the magnetar can be either a collapsar-type GRB/SN, 
a superluminous supernova (SLSN), or a neutron star merger event. All 
those cases are associated with significant optical emission. Moreover, if the 
progenitor is an off-axis GRB, the off-axis afterglow could last for weeks. 
Below we test if any of the expected optical emission is consistent with the 
upper limits obtained from the Pan-STARR data. 

For the case of SLSN, we choose SN2015bn as a representative and 
assume that it could serve as the progenitor of the magnetar in this study. 
SN 2015bn was located in a faint host galaxy at D; = 544.8 Mpc°!. The 
optical data are available at the Open Supernova Catalog®. We first fit the 
multi-wavelength optical light curves using the magnetar-powered model 
described in Ref. and obtain the physical parameters to describe the ob- 
servational data of SN 2015bn. As listed in Extended Data Table 4, after 
being updated on the luminosity distance and the magnetic field using those 
values of the hyper flare GRB 130310A, the parameter set is utilized to cal- 
culate the r-band magnitude using the same model. The resulted light curve 
is plotted as a solid green line in Figure 4. One can see the SLSN magnitude 
is way above the PAN-STARR upper limit. The weeks-scale coverage of 
PAN-STARR observations around To can rule out the existence of an SLSN, 

For the case of a collapsar-type GRB/SN, We choose SN 1998bw as 
an representative. SN 1998bw was located in a barred spiral galaxy (ESO 
184-G82) at Dz =~ 37 Mpc®. The data are also available in the Open Super- 
nova Catalog. After fitting the multi-wavelength optical light curves using 
the °°Ni-decay model™, the best-fit physical parameters expect for the lu- 
minosity distance (as listed in Extended Data Table 4) are used to calculate 
the r-band Type Ic SN light curve at 70 Mpc. The resulting light curve, as 
shown in Figure 4, is well above the PAN-STARR upper limit. We can thus 
also rule out the classical collapsar-type supernova origin of the magentar. 

For the case of a neutron star merger, one should consider the existence 
of a kilonova (KN). A comparable case is the observed AT2017gfo with 
GRB 170817A. AT2017gfo was located in an elliptical galaxy (NGC 4993) 
at Dz = 40 Mpc®. The optical data are available at the Open Kilonova Cat- 
alog”. We apply the model employed in Ref and fit the multi-wavelength 


‘https: //outerspace.stsci.edu/display/PANSTARRS/PS1+Stack+images 
Shttps://outerspace.stsci.edu/display/PANSTARRS/PS | + Image+Cutout+Service 
Shttps://sne.space 

Thttps://kilonova.space 


optical light curves to obtain the best-fit physical parameters as listed in Ex- 
tended Data Table 4. After updating the luminosity distance, we plot the 
r-band kilonova light curve at 70 Mpc in Figure 4. Our results suggest that 
the peak magnitude of the kilonova could be above the upper limit of point- 
source detection. On the other hand, the theoretical KN light curve fell in 
the observational gap of PAN-STARR. So the neutron star merger origin of 
the magnetar is allowed. 

Giving the fact that no GRB was found between ~ To — Tenax and To 
and that a kilonova could not be ruled out, the putative short GRB associated 
with the merger event was likely off-axis with respect to the observer, unless 
it was missed during the 35.4% off-FOV-time of Fermi/GBM. We compared 
the off-axis afterglow of GRB 170817A with the upper limit constrained 
from PAN-STARR observations. To do so, we first fit the r-band afterglow 
light curve of GRB 170817A using the structured jet model in Ref. 6, which 
gave us a set of best-fit physical parameters, as listed in Extended Data Ta- 
ble 4. After replacing the luminosity distance with 70 Mpc, we calculated 
the predicted off-axis afterglow in our case. Such a light curve is plotted 
in Figure 4. One can see that the off-axis afterglow of the short GRB is 
undetectable by PAN-STARR. We also calculated the on-axis short GRB af- 
terglow for the same set of parameters at 70 Mpc and found that such an 
afterglow can also survive the observational constraints if it occurred in the 
observational gaps. Considering the 35.4% off time for GRB, this suggests 
that an on-axis short GRB and afterglow as the magnetat progenitor is not 
ruled out. However, in view that the event rate density of off-axis neutron 
star merger events is much higher than that of the on-axis events, we regard 
an off-beam neutron star merger as a more likely progenitor of the magnetar. 


Late-time afterglow and multiwavelength counterparts 

The late-time afterglow of the event should undergo the slow-cooling 
phase. The flux density F, at frequency v follows F, œ y7?/?12-3)/4, Using 
the following typical parameters in a standard forward shock model: 


1. electron distribution power law index, p = 2.2 
2. energy fraction in electrons, € = 0.1 

3. energy fraction in magnetic fields, eg =0.01 
4. accelerated electron fraction, éy =1.0 


5. ISM density no = 1.0 cm? 


we can calculate the flux densities in X-ray at 10!8 Hz and optical band at 
10!4 Hz at t ~ 8.5 years as®: 


vey 10 ee 1 Pigs 2)/4 p(2+p)/4 (2-3p)/4 4-2 vg’ 
SAX) = 1.21 x 10 MY & 0168001 Eiso4o tgsyrs 427 , a1) 
and 
p-1 (p-2)/4.,(2+p)/4,(2-3p)/4 4-2 jaba 
f(O) = 3.04 x 10 my € 20.16B,0.01 Eiso49 'gsyrs 27% 14 ` (12) 


Magentasrs are also believed to be the energy source of Fast Radio 
Bursts (FRBs). For a typical FRB with luminosity of 10” erg/s, its ob- 
served flux density at 70 Mpc is about 50 Jy in typical raido bands, which 
is detected by the ground-based radio telescopes, such as CHIME, MeeKat, 
FaST and VLA. 

The luminosity of the X-ray bursts of the Galactic SGRs are typically 
~ 10° erg/s. Putting it at 70 Mpc, we found that it too faint for current 
X-ray missions including Chandra and Swift. 


Data Availability 

Processed data are presented in the tables and figures in the paper. 
Source and optical observational data are available upon reasonable re- 
quests to the corresponding authors. The Fermi/GBM data are publicly 
available at https: //heasarc.gsfc.nasa.gov/FTP/fermi/data. 


Code Availability 
Upon reasonable request, the code (mostly in Python) used to pro- 
duce the results and figures will be provided. 
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EXTENDED DATA 


Table 1 | Observational properties of the hyper flare event GRB 130310A. The total fluence and peak flux are calculated in 10-10,000 keV energy band. All errors 
correspond to the lo credible intervals. 


Observed Properties 


GRB 130310A 


Abrupt rise time 
Steep decay time 
Tog (burst) 

Too (precursor) 


~ 20 ms 

~ 900 ms 
2.90583 s 
1.307944 S 


Duration of the “sharp peak” ~ 27 ms 
Waiting time between precursor and main burst ~ 3.79s 

a at peak -1.15482% 

E, at peak 10.75443 MeV 
Time-integrated œ (main burst ) -1.10801 
Time-integrated E,(main burst) 2.731023 MeV 
Time-integrated œ (sharp peak ) —1.19 + 0.02 


Time-integrated E,(sharp peak) 
Total fluence 

Peak flux 

Possible host galaxy 

Luminosity distance 

Isotropic energy (precursor) Ey iso 
Isotropic energy (sharp peak) Ej iso 
Isotropic energy (main burst) Ey,iso 
Peak luminosity Ly 5 iso 

Average tail luminosity Lail 


9.80223 MeV 
+032 -5 2? 
Tock : a = os -1 
“U0 "0.44 ECES 

g0927191-170053 
69.41 Mpc (z=0.0155) 
+7.49 46 
5 61k 10 erg 
23718 x 10% = 
23002 x 105° ae -1 
390.25 gs 

5.63 x 10% erg s7! 


Spectral lag ~0s 
Minimal variability time scale ôt 0.0075 s 
f parameter 11.22 + 1.32 


Table 2 | Spectral properties of GRB 130310A 


Parameter(s) 
Component ti h Model a E, ot KT (keV) PGSTAT/dof 
Precursor -0.57 0.35 BB e 45.0677 307.2/349 
Sharp Peak 4.108 4.135 CPL -1.194007 9802.23#23888  354.3/349 
3.85 4.05 -2.094231 858.507 1908.71 163.3/349 
4.05 4.118 =1.29+0-08  3943,29+268162  270.2/349 
4.118 4.135 -1.15400 10749.48°7 57-2 356.3/349 
4.135 4.155 —1.04+003  4964.09+173392  306.5/349 
Main Burst 4.155 4.185 CPL 0.73190%  1153,87+?6238 275.7/349 
4.185 4.215 -0.790 p9 940.6171 29 251.1/349 
4.215 4.265 0.497010  593.03+1066 279.2/349 
4.265 4.485 0.69400 706.6373 60° 355.6/349 
4.485 5.25 -1.07199%  885.18+74726 316.8/349 
4.05 5.25 S110 70)  2732.60+232-15 421.2/349 
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Date Number of observations Observation time Exposure time(s) Filters Upper limit for a point source detection (mag) 
2012-11-04 1 14:46:23.563 30 y 18.72 
2012-11-25 2 15:28:02.527 30 Z Unconstrained 

15:43:20.493 30 Z Unconstrained 
2012-12-28 2 13:36:32.675 45 i 18.98 
13:54:38.657 45 i 18.94 
2012-12-29 2 14:09:08.155 40 T 19.49 
14:26:01.300 40 r 19.48 
2013-01-25 4 10:54:12.288 40 T Unconstrained 
10:55:56.784 40 r 19.38 
11:11:00.673 40 T Unconstrained 
11:12:44.948 40 r 19.38 
2013-01-26 2 11:30:33.887 45 i 19.06 
11:49:06.794 45 i 19.11 
2013-04-01 4 06:44:26.925 43 g Unconstrained 
06:51:45.444 43 g 20.87 
07:02:04.029 43 g Unconstrained 
07:09:21.358 43 g 20.75 
2013-04-24 5 05:32:38.923 30 y Unconstrained 
05:47:32.901 30 y Unconstrained 
05:54:35.440 30 y 18.91 
06:11:14.971 30 Z 18.97 
06:26:05.774 30 Z 18.91 
2013-11-18 1 15:52:03.043 80 y 18.45 


Table 3 | List of PAN-STARR observations that covered the host galaxy around 7). The up-limit calculation is detailed in Methods. 


SLSN Value Ic-SN Value Kilonova Value Off-axis afterglow Value 

Surface field Bs ~ 105 G Nickel fraction fyi = 0-023 Blue ejecta mass Mge =0.023 Mo fraction of magnetic energy eg = 10-4 
NS mass Mys = 1.79 Mo Opacity to y-rays ky = 0.11 cmg! Blue ejecta velocity vane, ~ 78000 km s~! fraction of electron population energy £e = 0.04 
Spin period Pepin = 1.60 ms Ejecta mass Mgj = 1.78 Mo Blue opacity blue =0.5cm2g7! power-law slope of the electron population p = 2.17 
Optical opacity «=0.20cem2g7! Host H number density nHhost ~ 10!7-8cm73 | Blue temperature blue — 3983 K ISM density ng = 0.002 cm~ 
Opacity to y-rays ky = 0.01 cm? g7! Photosphere temperature Tmin = 6607 K Red ejecta mass mgd = 0.050 Mo jet orientation 8y = 0.40 rad 
Ejecta mass Mej =0.1 Mọ Explosion time texp = -54.52 days Red ejecta velocity vga ~ 44700 km s7! jet core width 8e = 0.07 rad 
Ejecta velocity vej ~ 5460 km s7! Variance o = 0.38 mag Red temperature pred =3745K jet total width Oy = 0.47 rad 
Host galaxy Extinction Ay = 0.08 mag Ejecta velocity Vej ~ 77625 km s7! Red opacity eed = 10cm2g7! 

Temperature T=832x 103K Luminosity distance = 70 Mpc Variance o =0.24 mag 

Variance o = 0.18 mag Luminosity distance = 70 Mpe 

Luminosity distance = 70 Mpc 


Table 4 | Model parameters of the SLSN, Ic-SN, kilonova and off-axis afterglow. 
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Figure 1 | Light curve and duration of GRB 130310A. a, the black solid lines show the light curve obtained by the GBM-na data in energy range of 8-1000 KeV. The 
red solid line represents the level of background. The sharp spike (SP) is at ~ 4s. and marked with purple shaded area. The inset brackets the precursor region. b, black 
line represents the accumulated count light curve. The blue horizontal dashed (solid) lines are plotted at 5% (0%) and 95% (100%) of the total accumulated counts, 
respectively. The two regions marked by the green vertical dashed lines are the Too intervals of the precursor and the main burst. 
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Figure 2 | Multi-wavelength light curves, high-energy photon detection and time . a, Multi-wavelength light curves detected with GBM and LAT. The filled blue 
circle is an event detected by LAT around the peak time of the GBM light curve. b, LAT events with energy above 100 MeV in 0 — 800 s. Different colors represent 
different ranges of probabilities that the photon count may originate from this event. c, Comparison between LAT-LLE and GBM light curves. d, Energy dependent 
lags between the lowest energy(10-20 keV) band and any higher energy bands. All error bars represent 10 uncertainties. 
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Figure 3 | Spectral evolution of the GRB 130310A. The time-resolved spectra are fitted with a CPL model for the main burst and a black body model for the 
precursor. The top panel shows the evolution of the low-energy photon index, œ for the main burst. The bottom panel shows the evolution spectral peak, Ep, for both 
the precursor and the main burst, where an equivalent Ep ~ 2.82 kT (Ref. ©) is adapted for the precursor. All error bars represent 1-c uncertainties 
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Figure 4 | Too — f(fer) plot. The blue circles represent the f value of short GRBs, and black circles represent the fe of long GRBs. The effective amplitude parameter 
ter = ~ for a long GRB is measured by scaling down the burst until “pseudo GRB” is shorter than 2 s (Ref.**) .The green vertical dashed line is the boundary of 2 s. 
The green diamond represents the MGF GRB 200415A. The red-star-marker highlights the sharp peak of GRB 130310A . 
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Figure 5 | Localization Test Statistics (LTS) map of the hyper flare GRB 130310A obtained with LAT data in the time range of 0-800 s. The white, orange, 
green and black counters are the 68.3%, 99.5%, 99.7% and 99.9% confidence regions, respectively. 30 IPN error box !? (blue parallelogram) and the potential host 
galaxy g0927191-170053 (red star) are over-plotted on top of the LAT error circles. 
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Figure 6 | An example illustrating the method of the up-limit calculation. a, A warp image of the host galaxy, g0927191-170053, obtained from the PAN-STARR 
archive for the observation performed on 2012-12-28T13:36:32.675. The green ellipse annuli is the region to measure the up-limit of the point-source transient. b, The 
magnitude distribution of each pixel within the annuli. The black solid line is the probability density function, and the black dashed line marks the 3 level. 
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Figure 7 | The wavelet plot of the detrended data. A significant signal is shown at 12 HZ between 4.4 and 5.0 s. The left plot is derived from WWZ (Weighted 
Wavelet Z-Transform) transform of the detrended data. The right panel shows the mean wwz power for different frequencies. 
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Figure 8 | An illustration shows that the periodic signal can become undetectable depending on its modulation strength to the observed flux. In each plot, the 
simulated light curve is composed of a FRED-shaped pulse (solid red line on the top panel), a random noise signal (blue in the middle panel), and a 12.4-Hz periodic 
signal applied between 4.4 s and 7.0 s (solid red line in middle panel). With the strength of the periodic signal decreasing from a to f (see Methods), one can see the 
detectable range of the 12-Hz period becomes more concentrated to its maximum-strength region. The detecting range of the periodic signal in d is similar to that of 


GRB 130310A. 
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